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ABSTRACT 

We use two 4.5/im Spitzer (IRAC) maps of the NGC 1333 region taken over a ~ 7 yr interval to 
determine proper motions of its associated outflows. This is a first, successful attempt at obtaining 
proper motions of stellars outflow from Spitzer observations. For the outflow formed by the Herbig- 
Haro objects HH7, 8 and 10, we find proper motions of ~ 9-13 km s~^, which are consistent with 
previously determined optical proper motions of these objects. We We determine proper motions for 
a total of 8 outflows, ranging from ~ 10 to 100 km s~^. The derived proper motions show that out of 
these 8 outflows, 3 have tangential velocities < 20 km s~^. This result shows that a large fraction of 
the observed outflows have low intrinsic velocities, and that the low proper motions are not merely a 
projection effect. 

Subject headings: ISM: jets and outflows — ISM: Herbig-Haro objects — ISM: star formation — cir- 
cumstellar matter — stars: formation — infrared: ISM — ISM: individual objects 
(HH7-11) — ISM: individual objects (NGC 1333) 



1. INTRODUCTION 

Starting with the work of Herbig & Jones (1981), who 
studied the HH 1/2 outflow, the study of proper mo- 
tions of Herbig-Haro (HH) objects has become one of 
the fundamental tools in exploring the properties of out- 
flows from young stars. Proper motions of young stel- 
lar outflows were first measured at radio wavelengths by 
Rodriguez et al. (1989), who studied the Serpens "triple 
source". Later, Noriega- Crespo et al. (1997) first ob- 
tained proper motions of outflows at IR wavelengths of 
HHl, from H2 2.12/im images. 

In contrast to what has happened at optical and radio 
wavelengths, the study of IR proper motions of young 
stellar outflows has not flourished. This is probably due 
to the fact that older IR images typically cover only rel- 
atively small fields, so that it is many times impossible 
to obtain a good centering/scaling between them and 
newer (larger field) IR images. Previous IR proper mo- 
tion studies include Noriega- Crespo et al. (1997: HHl), 
Chrysostomou et al. (2000: HH7-11, 25-26 and 33/40) 
and Raines et al. (2000: GGD37). 

In this paper we make a first attempt to measure 
proper motions from Spitzer-IRAC IR images of young 
stellar outflows. This is possible because even though 
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the IRAC pixels are relatively large (l'^2), some outflows 
have now been observed over time intervals of ~ 10 yr. 
For example, for an outflow with plane-of-the-sky veloci- 
ties of ~ 100 km s~^ at a distance of ~ 500 pc, offsets of 
~ 0.35 pixels are expected over ~ 10 yr. As shown ini- 
tially by Raga, Barnes & Mateo (1990), such offsets can 
be measured with reasonable accuracy in CCD images of 
stellar outflows. 

There are three fundamental reasons that allow us to 
use IRAC to measure proper motions in the mid-IR, de- 
spite its relative low angular resolution (FWHM~2'') in 
comparison with optical or near-IR ground based obser- 
vations. First of all is the outstanding stability of the 
pointing control system of the Spitzer Space Telescope 
over the past nine years. Secondly, the unchanged per- 
formance of its optical system and its detectors. Thirdly, 
the ability of IRAC to map large regions of the sky, thus 
providing plenty of point sources that can be used to 
create and fix a reliable reference frame to compare ob- 
servations from different epochs. 

We study two 4.5/im IRAC frames of the NGC 1333 
region, obtained in February /September 2004 and Octo- 
ber/November 2011. These observations are described in 
detail in section 2. 

These images allow us to obtain the proper motions of 
some of the HH7-11 knots (HH7, 8 and 10), which can 
be compared with previous optical (Herbig & Jones 1983; 
Noriega- Crespo & Garnavich 2001) and IR (Chrysosto- 



2 



Raga, Noriega- Crespo, Carey & Arce 



mou et al. 2000) proper raotion Hieasureraents of this 
outflow. This is done in section 3. Also, these images al- 
low us to obtain the proper motions of other outflows in 
NGC 1333, as described in section 4. Finally, our results 
are summarized in section 5. 

2. OBSERVATIONS 

NGC 1333 was one of the first star forming regions ob- 
served with the Spitzer Space Telescope (Werner et al. 
2004) with the Infrared Array Camera (IRAC) (Fazio et 
al. 2004) during the Cryo-mission as part of the IRAC 
GTO program of its principal investigator (PI) (PID 6, 
G. Fazio). This first epoch was complemented approx- 
imately six months later by the Cores-to-Disks Spitzer 
Legacy program as part of their mapping of the Perseus 
molecular cloud (PID 178, PI N. Evans). 

Multiple epoch observations of star forming regions rel- 
atively near the ecliptical plane are essential to weed out 
asteroids when trying to identify real proto-stellar ob- 
jects (see e.g. Evans et al. 2009). The more recent 
data were obtained in October-November 2011, and are 
being currently analyzed to study the time variability 
of young stellar objects (YSOs) in the region, as part 
of the Warm-mission Spitzer Exploration program YSO- 
VAR (PID 61026, PI J. Stauffer). YSOVAR has observed 
the core of NGC 1333, a ~ 16'xll' region centered on 
03h29ml0.175s -h31dl6ml6.0s, over a period of nearly a 
month using 73 different visits. Each visit lasted approx- 
imately 15 minutes with a mean coverage of ~ 4 — 5 in 
the final maps, although some pixels in the overlap re- 
gions can double these values. Our final 2010 image has 
a mean coverage of ~ 300, i.e., on average each pixel has 
seen the object some 300 times, and since the observing 
time per pixel is 12 sec, this corresponds to a ~ one hour 
total exposure time. 

This central map has determined the size of the region 
that we have used for our analysis of the proper motions 
between the two (2004 & 2011) available epochs. The 
final map is actually a bit smaller to avoid the scattered 
light present around the very bright young star LZK 12 in 
the very deep YSOVAR image. The Cryo-mission obser- 
vations were carried using all IRAC channels (3.6, 4.5, 
5.8 & 8.0/im), while in the Warm-mission observations 
we can only use channels 1 and 2 (3.6 and 4.5/im). For 
the three mentioned programs, the data was collected us- 
ing the High-Dynamic-Range (HDR) mode with a 12sec 
integration time for the 'long' frames and 0.6sec for the 
'short' ones. For the purpose of our analysis we have con- 
centrated on the data at 4.5/im, which for young stellar 
outflows has proven to be a very good tracer of the molec- 
ular Hydrogen emission of some of the bright pure rota- 
tional transitions [(0-0 S(9) at 4.6947, 0-0 S(10) at 4.4096 
and 0-0 S(ll) at 4.1810/im], and therefore, of their shock 
excited emission (Noriega-Crespo et al. 2004a,b; Looney 
et al. 2007; Tobin et al. 2007; Ybarra & Lada 2009; 
Maret et al. 2009; Raga et al. 2011; Noriega-Crespo & 
Raga 2012). 

For the IRAC observations we have used the latest ver- 
sions of the data, S18.18.0 (Cryo) and S19.0.0 (Warm), 
and a scale of 0.3'Ypix for the final maps, i. e. a seventh 
of the IRAC FWHM-- 2" angular resolution. The maps 
of the different epochs are created with the same fiducial 
frame, using MOPEX (Makovoz, Khan & Masci 2006), 
which brings them very close to the same frame of refer- 




Figure 1. Top frame: 2004 — 2011 difference image. Center: 2011 
image. Bottom: 2004 image. In the difference image, the knots 
along outflows with proper motions show up as "positive/negative" 
double ridge structures, and the white disks are masked stars. 

ence. A dozen non-saturated stars common between the 
two epochs are then used for the final alignment of the 
images, with a typical centroid accuracy of less than a 
tenth of a pixel. The two resulting maps (2004 and 2011) 
are shown in Figure 1. 

Figure 1 shows the two final images (2004 and 2011), 
and the difference between these two images. In the 
difference map, one can clearly see that some of the 
compact emission features show clear proper motions. 
For example, a bow-shaped region in the SW region of 
the frame shows a "positive/negative" ridge combination 
which clearly indicates a motion of this feature. Such 
motions are analyzed in detail below. 

Although there are no major changes in the optics of 
the telescope or the performance of the detectors be- 
tween the IRAC Cryo and Warm missions, a comparison 
between the 2004 and 2011 images shows stronger opti- 
cal ghosts around the bright sources in the Warm data. 
This effect is probably enhanced because of the very deep 
coverage of the YSOVAR data (with ~ 300 exposures 
of 12 sec for each pixel, see above). Figure 2 shows a 
3.5'x3.5'region of NGC 1333 with some bright sources, 
where one can see the typical 'heart shape' IRAC Point 
Spread Function, with the nearby "triangular ghosts" 
marked with the white boxes (12^^x12''). These are an 
annoyance more than a problem in measuring the proper 
motions, since they appear as real time variations be- 
tween epochs. For this reason we have masked some of 
them in both final (2004 and 2011) images. The obser- 
vations are summarized in Table 1. 

3. THE HH7-11 OUTFLOW 

In Figure 3, we show the 2004 frame of the HH7-11 
outflow. HH 9 is not seen in the 4.5/im IRAC maps, and 
HHU is lost within the PSF of SVS13. The image shows 
a knot 2^" to the S of HHIO, which we have labeled 
"knot A" (see Figure 3). 

Figure 3 also shows proper motion vectors calculated 
with the "proper motion mapping" technique described 
by Raga et al. (2012). This (rather simplistic) technique 
is described in Appendix A. In particular, the HH7-11 
proper motions shown in Figure 3 were computed with 
cross correlations in boxes of size L = 30 pixels (with 
each pixel corresponding to 0''.3, see section 2) which 
have peak fluxes of at least frnin = 2.0 mJy/sterad (see 
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Figure 4. 2004 frame (displayed with the logarithmic scale given 
in mJy/sterad by the top bar) of the studied field. The dark, circu- 
lar regions correspond to masked stars. The blue arrows indicate 
the axes of the outfiows described by Davis et al. (2008), which we 
have numbered 1 through 9. The yellow boxes (labeled A through 
F) indicate interesting subfields, which are shown in an expanded 
scale in Figures 5-10. The white arrows indicate the proper motion 
velocities (for a distance of 220 pc) obtained through our "cross 
correlation mapping" technique. N is up and E to the left. 



Figure 3. The map shows HH7-11 outfiow lobe in the 2004 map. 
The white arrows indicate the proper motion velocities (for a dis- 
tance of 220 pc to the system) computed within the cross corre- 
lation boxes with high enough intensities. We detect motions for 
HH 7, 8 and 10 (HH 11 being confused with the stellar PFS off the 
NW corner of the frame). N is up and E to the left. The image 
is displayed with a logarithmic scale given (in mJy/sterad) by the 
top bar. 

Appendix A). The proper motion velocities were calcu- 
lated assuming a distance of 220 pc to NGC 1333 (Cernis 
1990; Hirota et al. 2008). 

The HH7, 8 and 10 knots have velocities ^ (10 ± 
4) km s-i (HH7: 13 km s-^ HH8: 10 km s-^ HHIO: 9 
km s~^), with HH7 and HH8 moving approximately par- 
allel to the outflow axis, and HHIO showing a possibly 
significant deviation towards the S. These proper motions 
were obtained by averaging the proper motions from the 
two boxes closest to the peak emisison of the knots. Also 
the region between HH7 and 8 shows filamentary struc- 
tures which apparently share the kinematics of these two 
HH knots (see Figure 3). 

Knot A has a significantly higher proper motion veloc- 
ity of (25 ± 4) km s~^, pointing approximately E. This 
knot could belong to one of the other outflows present in 
this complex region (see section 4 and Davis et al. 2011). 

Finally, relatively large proper motions are seen within 



~ 20^' from the SVS13 outflow source (outside the upper, 
right hand side corner of Figure 3), which are probably 
related with the spikes of the PSF, and are unlikely to 
represent real motions in the flow. 

The proper motions that we have obtained for HH 7, 
8 and 10 are roughly consistent with the velocities of 
~ (20 ± 10) km s~^ determined for these knots with the 
proper motion determinations of Noriega- Crespo & Gar- 
navich (2001), scaled to a distance of 220 pc to HH7-11. 
The proper motions of Herbig & Jones (1981) are also 
consistent with these relatively low velocities. The ve- 
locities of ~ 400 km s~^ determined by Chrysostomou et 
al. (2000) are inconsistent with the three other indepen- 
dent (optical and IR) proper motion determinations. 

4. OTHER OUTFLOWS 

In Figure 4, we show a 450 x 550^' field (covered in 
the two epochs) approximately centered on the HH 7-11 
outflow. On the 2004 frame, we show the axes of the 
molecular outflows (which we have numbered 1 through 
9) described by Davis et al. (2008). As already seen in 
the diffenerce image between the two epochs (see Figure 
1 and section 2), there are several regions of extended 
emission which show detectable proper motions. 

The proper motions shown in Figure 4 correspond to 
cross correlations within L — 60 (18^^) pixel boxes of the 
two epochs, with fluxes of at least fmin = 0.9 mJy/sterad 
(see Appendix A). 

We have chosen 6 regions (labeled A through F in Fig- 
ure 4) which show the more interesting proper motion 
detections. These regions are described in the following 
subsections. 
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Figure 5. Subfield A of Figure 4. The 2004 map is shown together 
with the proper motion velocities. N is up and E to the left. 

4.1. Region A 

This region (see Figure 5) includes the HH 7-11 outflow 
(the smaller box within region A in Figure 4, correspond- 
ing to the blue-shifted, 7-11 lobe, is shown in Figure 3). 
It is clear that while organized proper motions (aligned 
with the outflow axis) are seen along the HH 7, 8, 10 
knots, such motions are not seen in the NW outflow lobe. 

To the S of the SVS 13 source, relatively large proper 
motions with Easterly direction are seen, but it is not 
clear that these motions are associated with well defined 
knots (except for "knot A" of Figure 3, see section 3). 
These motions might be associated with a previously un- 
detected outflow. 

In the center of the frame, at ~ 30^' from the Southern 
edge of the frame, we see a bow-like region with North- 
ward proper motions of ~ 10 km s~^. These marginally 
detected motions could be associated with the Northern 
lobe of outflow 4 (see Figure 4). 

Close to the SW corner of the frame, we see a region 
with a broken-up, bow-like structure with proper motion 
velocities of ~ 20-70 km s~^ to the N or NE. This motion 
appears to be shared by some of the knots seen in the 
NW region of Figure 5. All of these Northward directed 
knots appear to belong to outflow 6 (see Figure 4). 

4.2. Region B 

This region (see Figures 4 and 6) shows two systems 
of knots with different proper motions. One set of knots 
(in the Northern part of Figure 6) has proper motion 
velocities of 20-100 km s~^ directed approximately to 
the SSE. This set of knots appears to be associated with 
outflow 2 (see Figure 4). 

A second set of knots (in the Southern half of Figure 6) 
shows proper motions of similar magnitudes (10-60 km 
s~^), but directed Eastwards. These knots could be as- 
sociated with outflow 8 (see Figure 4), though their mo- 
tion would imply a substantial deflection from the out- 
flow axis. Alternatively, these knots could belong to a 
previously undetected outflow. 

4.3. Region C 

This region (see Figures 4 and 7) shows knots with 
proper motions of 10-60 km s~^ directed towards the 
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Figure 6. Subfield B of Figure 4. The 2004 map is shown together 
with the proper motion velocities. N is up and E to the left. 

SE. These knots are therefore likely to be a prolongation 
of the S lobe of outflow 2 (see Figure 4). 

4.4. Region D 

The Southern part of this region (Figure 8) shows a 
bow-like structure with Eastwards proper motions of 20- 
40 km s~^, apparently associated with the E lobe of out- 
flow 3 (see Figure 4). In the NW corner of the frame, 
we find a set of knots with Southward proper motions 
of 20-40 km s~^, which are likely to be associated with 
outflow 5 (see Figure 4). 

The region which is morphologically aligned with out- 
flow 8 (extending from the N center to the SW corner of 
Figure 8) does not show proper motions aligned with the 
outflow axis. 

4.5. Region E 

In the SW corner of this region (Figure 9) we see a set 
of knots with Western proper motions of 10-80 km s~^ 
which are morphologically and kinematically associated 
with the W lobe of outflow 3 (see Figure 4). 

In the SE corner, we see a bow shaped structure with 
Southern proper motions of up to 100 km s~^, which is 
likely to be associated with outflow 4. In the center of 
this region we see knots with Southern proper motions 
of ~ 60 km s~^ which might be associated with outflow 
6 (see Figure 4). 

Finally, in the NE corner of this field (see Figure 9) 
we see aligned proper motion vectors towards the E, cor- 
responding to velocities ~ 15 km s~^. These motions 
could be associated with the Eastern lobe of outflow 7 
(see Figure 4). 

4.6. Region F 

This region (Figure 10) shows some high Southwards 
proper motions of ~ 50 km s~^ (actually, one of the 
Southwards proper motion velocities is > 100 km s~^, 
but for a cross-correlation box without any clear feature), 
which do not have a clear identification with the outflows 
of Davis et al. (2008). 

Also, this region might be showing evidence of a bipo- 
lar structure (with oppositely directed proper motion ve- 
locities of ~ 10 km s~^, indicated with two open circles 
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Figure 7. Subfield C of Figure 4. The 2004 map is shown together 
with the proper motion velocities. N is up and E to the left. 
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Figure 8. Subfield D of Figure 4. The 2004 map is shown together 
with the proper motion velocities. N is up and E to the left. 

in Figure 10), which might be associated with outflow 9 
(see Figure 4). 

4.7. Proper motions for the outflows of Davis et al. 

(2008) 

As described above, we detect proper motions for knots 
which are hkely to be associated with several of the out- 
flows of Davis et al. (2008). These outflows are numbered 
1 through 9 in Figure 4. 

For each of these outflows, the following results are 
obtained: 

• outflow 1 - we measure proper motions of up to 20 
km s~^ along the SE lobe of this outflow. These 
proper motions are consistent with the optical 
proper motions measured for HH 7-11 by Noriega- 
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Figure 9. Subfield E of Figure 4. The 2004 map is shown together 
with the proper motion velocities. N is up and E to the left. 

Crespo & Garnavich (2001). We do not see orga- 
nized motions in the region that would correspond 
to the NW lobe of this outflow, 

• outflow 2 - we detect two sets of knots (within re- 
gions B and C of Figure 4, respectively) with proper 
motions of up to ~ 100 km s~^ aligned with the SE 
lobe of this outflow, 

• outflow 3 - we detect proper motions along both 
the E (region D) and W (region E) lobes of this 
outflow, with velocities of up to ~ 80 km s~^, 

• outflow ^ - we detect a Southern bow shock with 
proper motions of ~ 100 km s~^ (region E) and 
possibly a Northern bow shock with much lower 
(~ 10 km s~^) velocities (region A), which might 
be associated with this outflow, 

• outflow 5 - we detect knots with velocities of ~ 30 
km s~^ which might be associated with the South- 
ern lobe of this flow (region D), 

• outflow 6 - we detect knots with proper motions of 
~ 50 km s~^ which are likely to be associated with 
the Northern (region A) and Southern (region E) 
lobes of this outflow, 

• outflow 7 - the E lobe of this outflow coincides with 
IR knots with velocities of ~ 15 km s~^, 

• outflow 8 - we do not obtain proper motion veloc- 
ities aligned with this outflow, 

• outflow 9 - we marginally detect proper motion ve- 
locities of ~ 10 km s~^ which might be associated 
with this outflow. 

Therefore, we detect knots with measurable proper mo- 
tions which are likely to be associated with all but one 
of the outflows of Davis et al. (2008). 
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Figure 10. Subfield F of Figure 4. The 2004 map is shown to- 
gether with the proper motion velocities. The two open circles 
indicate the regions that could be associated with the two lobes of 
outflow 9 (see section 4.6). N is up and E to the left. 

4.8. The proper motion velocity distribution function 

If we are detecting most of the outflows in NGC 1333, 
the measured proper motion velocities vt allow us to 
obtain the distribution function associated with the out- 
flows from a group of stars currently being formed. If we 
consider the maximum value of vt for the 8 outflows of 
Davis et al. (2008) with measured proper motions (see 
subsection 4.7), we find that we have 3 with maximum 
Vt in the ^ 20 km s~^ range, 1 outflow in each of the 
20 40, 40 60 and 60 ^ 80 km s-^ ranges, and 2 
outflows with Vt in the 80 100 km s~^ range. The 
resulting distribution function (normalized to a total of 
5 outflows, see below) is shown in Figure 11. 

The distribution function shows an apparent peak at 
low proper motion velocities (vt), and possibly also a 
marginal peak in the higher velocity bin. One can show 
that this distribution is not consistent with a set of out- 
flows with identical intrinsic velocities Vj and with arbi- 
trary orientation angles with respect to the plane of the 
sky. 

It is straightforward to show that if we have a set 
of outflows with outflow velocity Vj and with arbitrary 
orientations, the expected distribution function for the 
proper motion velocities vt is: 



f{VT) = 



1 



-1/2 



(1) 



This distribution function is plotted in Figure 11, consid- 
ering Vj = 100 km s~^. In this Figure, the distribution 
function obtained from the 8 observed outflows has been 
normalized to 5 outflows, i.e., not considering the 3 out- 
flows in the ^ 20 km s~^ bin. A comparison between 
the observed and theoretical distributions shows that a 
single population of vj = 100 km s~^ outflows (at random 
orientations) does not produce a proper motion velocity 
distribution consistent with the data. 

Particularly, from Figure 11 it is clear that the ob- 
served outflows have a large population in the lowest ve- 
locity bin (0 20 km s~^), for which the ^^Vj = 100 km 
s~^ model" predicts low frequencies. This large popula- 
tion of the lowest proper motion velocity bin therefore 
implies that some of the outflows in our group indeed 




have a low, full spatial velocity, of at most Vj 

-1 



20 km 



[km s ^] 

Figure 11. The dashed histogram shows the distribution function 
obtained by binning the peak proper motion velocities of the 8 
outflows of Davis et al. (2008) which we have detected, normalized 
to a total of only 5 outflows. The solid curve is the distribution 
function expected for a set of outflows of intrinsic velocity vj = 
100 km s"-*^ at random orientations (see equation 1). 

5. SUMMARY AND CONCLUSIONS 

In this paper, we present a first attempt to derive 
proper motions of young stellar outflows from Spitzer 
(IRAC) images. Our analysis of the NGC 1333 region 
(which includes the HH 7-11 outflow) shows that it is 
indeed possible to do this, which opens up the possibil- 
ity of obtaining new "second epoch" observations with 
Spitzer, in order to derive proper motions for a relatively 
larger number of IR outflows. 

We have analyzed two epochs at a time interval of ^ 7 
yr, and the fact that proper motions are indeed detected 
can be seen by taking the difference between the two 
frames (see Figure 1). Using these two frames, we first 
derive the proper motions of HH 7, 8 and 10. 

The resulting proper motions correspond to plane of 
the sky velocities of ^ (10 =b 4) km s~-^ (see Figure 3 
and section 3), at a distance of 220 pc, which are con- 
sistent with the optical proper motions of these objects 
derived by Herbig & Jones (1983) and Noriega-Crespo 
& Garnavich (2001). Our results do not agree with the 
~ 400 km s~^ IR proper motions of Chrysostomou et al. 
(2000). 

We have then derived proper motion velocities for the 
emission knots in a 450 x 550^' field, many of which 
are aligned with the outflows described by Davis et al. 
(2008). For one of these outflows (outflow 3, see Figure 
4), we find a clear bipolar motion (with velocities of ~ 80 
km s~^). The CO emission of this outflow is clearly seen 
in the observations of Curtis et al. (2010). Another re- 
markable result is the ~ 100 km s~^ velocity derived for 
the Southern bow shock of outflow 3 (see Figures 4 and 
9). 

From the proper motions of the 8 outflows of Davis et 
al. (2008) for which we have measurements, we obtain 
a proper motion velocity distribution function (see sec- 
tion 4.8). We find that this distribution function is not 
consistent with a set of outflows with identical outflow 
velocities (of ~ 100 km s~^) at random orientations. The 
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empirically derived distribution function indicates that 
some of the observed outflows have intrinsic outflow ve- 
locities Vj 100 km s~^, and that others have a much 
lower, Vj ~ 10-20 km s~-^ velocity. This consists of a 
first suggestion that the velocities vj (of a group of out- 
flows from recently formed stars in the same cloud) have 
a bimodal distribution. 

Even though the derived proper motions only give par- 
tial constraints on the kinematics of the outflows in NGC 
1333, they do provide some of the information necessary 
for modelling this system. This is of particular interest, 
since it appears to be a case in which the outflows have a 
strong influence on the structure of their parental cloud, 
and might even be triggering the formation of more stars 
(see Sandell & Knee 2001). Some other clouds might also 
be in a regime of strong perturbation from the outflows 
from YSOs (see, e.g.. Are et al. 2010; Duarte-Cabral et 
al. 2012). 

It would be interesting to use the kinematics derived 
from our proper motion determinations to compute simu- 
lations of jet-driven turbulence such as the ones of Naka- 
mura & Li (2007) and CarroU et al. (2010). Would we 
then obtain a perturbed cloud such as NGC 1333 (Sandell 
& Knee 2001; Arce et al. 2010)? This would be an inter- 
esting complement to the quite extensive work that has 
been done on this object (see the papers quoted above 
and Quillen et al 2005; Maret et al. 2009; Padoan et al. 
2009; Arce et al. 2010). 
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APPENDIX 

Appendix A: Proper motion mapping 

In order to derive proper motions from pairs of CCD images it has now become standard to define boxes including 
emitting "knots", and to carry out cross-correlations (between the 2 frames) of the emission within the boxes. The 
proper motion is then obtained from a fit to the peak of the cross-correlation function. This method has proven to be 
better than carrying out direct fits (of e.g., a Gaussian or a 2D paraboloid) to the observed emission features, because 
the cross-correlation functions (which are integrals of the emission within the chosen boxes) have higher signal-to- noise 
ratios than the images. 

For images with many knots, it is simpler to use a regular array of boxes (or "tiles", as called by Szyszka et al. 2011) 
for computing cross correlations. Proper motions for the emission within each of the boxes are then obtained from fits 
to their respective cross-correlation functions. 

In this paper, we use tha setup described by Raga et al. (2012). We define square boxes of side L, with spacings 
of I//2 along each of the coordinates of the images. Therefore, neighbouring boxes have a superposition region of size 
L/2 along each axis, as shown in the schematic diagram of Figure 12. 

For a given box of size L(see Figure 12), we first check whether or not the condition / > fmin (where / is the 
intensity of the image and fmin is a user specified "minimum fiux") is satisfied in at least one pixel within a central 
"inner box" of size L/2. If this condition is met at least for a single pixel in each of the two epochs that are being 
analyzed, the cross-correlation function (within the L-size box in the two frames) and the proper motion (from a 
paraboloidal fit to 3 x 3 pixels centered on the peak of the cross correlation function) are computed. 

Therefore, in order to compute a "cross correlation map" , it is then only necessary to specify the box size "L" and 
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Figure 12. Schematic diagram showing one of the cross correlation boxes (thick hne box, of size L). The sohd circle is the central 
position of this box, and the open circles are the central positions of the nearby cross correlation boxes. The inner box (drawn with a thin 
line) of size L/2 is used for evaluating whether or not the peak flux satisfies the criterion for calculating a proper motion 



Table 1 

NGC 1333 Observations 



Program ID Request Key Observation Time PI 

6 3652864 2004-02-10 Fazio, G. 

178 5793280 2004-09-08 Evans, N. 

61026^ 29323008 2011-10-12 Stauffer, J. 

29319424 2011-11-13 



^ There are 73 independent observations (Request Key) in this program; the first and last are shown. 

the minimum flux fmin- Actually, there is another "hidden" parameter which is the centering of the cuadricule of 
boxes with respect to the pixels of the image. This is fixed arbitrarily so that the first pixel of the frames that are 
being studied coincides with the edge of one of the cross- correlation boxes. 

We find that the peaks of the cross correlation boxes can be fitted to within ~ 0.1 pixel, a precision similar to the 
one of the alignment between the two epochs that we have used. Therefore, the formal errors of our proper motion 
determinations are of 0.1 pixel, corresponding to a velocity of ~ 4 km s~^. 



